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O . Abstract 

As a possible source of the single transverse spin asymmetry, we study the contribution from 
^ ■ purely gluonic correlation represented by the twist-3 "three-gluon correlation" functions in 

j^ . the transversely polarized nucleon. We first define a complete set of the relevant three-gluon 

correlation functions, and then derive its contribution to the twist-3 single-spin-dependent 
cross section for the D-meson production in semi-inclusive deep inelastic scattering, which 
is relevant to determine the three-gluon correlations. Our cross-section formula differs from 
the corresponding result in the literature, and the origin of the discrepancy is clarified. 



1 Introduction 

Clarifying the origin of the large single spin asymmetries (SSAs) observed in various high- 
energy semi-inclusive processes [l]-[7] has been a big challenge during the past decades. (See 
[8] for a review.) They are now understood as direct consequences of the orbital motion of 
quarks and gluons and/or the multi-parton correlations inside the hadrons, and thus provide 
a new opportunity for revealing the QCD dynamics and hadron structures, which do not 
appear in the parton models and perturbative QCD at the conventional twist-2 level. Satis- 
factory formulation of these effects in QCD, however, requires sophisticated framework for 
which lots of technical developments are involved. Up to now, the SSAs have been formu- 
lated based on the (naively) "T-odd" distribution/fragmentation functions [9]- [18] within 
the transverse-momentum-dependent (TMD) factorization approach [19, 20, 21], and on the 
twist-3 multi-parton correlation functions in the collinear factorization approach [22]- [38]. 
These two mechanisms, in principle, cover different kinematic regions, such that the former 
approach is suitable for treating SSAs in the region of the small transverse momentum of a 
particle observed in the final-state, while the latter is designed for systematic description of 
SSAs at large values of the corresponding transverse momentum. On the other hand, in the 
intermediate region of the corresponding transverse momentum where both approaches are 
valid, it's been shown for some structure functions in semi-inclusive deep inelastic scattering 
(SIDIS) and in the Drell-Yan cross section that the two approaches provide an equivalent 
description of SSAs [39, 40, 41, 42]. Therefore, in practice, the two approaches, tied together 
at the intermediate transverse-momentum region, can be regarded as providing a unique 
QCD effect leading to SSAs over the entire kinematic regions. Although our understanding 
on SSAs has made a great progress armed with these mechanisms, further studies are yet 
to be done for a complete clarification of all the effects responsible for SSAs. Among such 
effects, in particular, the role of purely gluonic effects has not been widely studied in the 
literature. Since the gluons are ample in the nucleon, they are potentially an important 
source of SSAs. 

In this paper, we study a purely gluonic effect as an origin of SSAs in the framework of 
the collinear factorization. To this end, we work on SSAs for a heavy meson production in 
SIDIS, in particular, the Z)-meson production, ep — )■ eDX. Since this process is induced 
dominantly by the cc-pair creation through photon-gluon fusion, this is the most relevant 
process to probe the gluonic effects for SSAs, together with D-meson production in pp 
collisions [43, 33] ongoing at RHIC [44]. In the collinear factorization framework, SSA is a 
twist-3 observable and thus the gluonic effects responsible for SSAs are represented by the 
twist-3 gluon correlation functions in the polarized nucleon, which were first introduced in 
the most general form by Ji[24]. Recently the authors of [32] studied the SSAs in SIDIS, 
ep^ — 7- eDX, applying the twist-3 mechanism. They derived a formula, in the leading 
order with respect to the QCD coupling constant, for the contribution to the single-spin- 
dependent cross section from a "three-gluon" correlation function. In our opinion, however, 
the starting formula used to derive the cross section in [32] is incomplete and cannot lead 
to the correct result. So we revisit the same issue in this paper. 

As is well-known, a twist-2 gluon distribution inside the nucleon is defined in terms 



of the gauge-invariant lightcone correlation function of the gluon field-strength tensors, 
schematically written as {F""^ F^^) . Likewise, the twist-3 "three-gluon distribution" func- 
tions are defined through the gauge-invariant correlation functions {C^'^F^^FJ^ F2^)i "with 
the structure constants of the color SU(3) group, C"^'^ = if"-'"^ and C'^'^ = d"''"^. On the other 
hand, when we derive the three-gluon contribution to the cross section for ep^ — )■ eDX with 
the large transverse momentum of the D meson, we start our analysis with a cut forward- 
amplitude for the cross section, in which the correlation functions of the gluon fields appear 
in the form ~ (A^A^AJ). Extraction of the twist-3 effect relevant for SSAs from the cor- 
responding amplitude, converting eventually the associated nucleon matrix elements from 
(A^A^AJ) into the gauge-invariant forms {C^^ F^^ FJ^^ F^^) ^ is a highly nontrivial issue, un- 
like straightforward calculations of the twist-2 cross sections. In this connection, it is worth 
mentioning that, for the case where the twist-3 quark-gluon correlation functions partici- 
pate, the necessary formulation was achieved in [27]; there. Ward identities played a crucial 
role to prove the factorization property and gauge invariance of the corresponding twist- 
3 cross section. Similarly, we will show that, owing to the constraints from the tree-level 
Ward identities satisfied by the relevant partonic hard parts, the twist-3 contribution to the 
cross section, associated with the three-gluon correlation functions {A^AIIA].)^ can be recast 
into the factorized expression in terms of the gauge-invariant functions {C^'^F^^FJ^ F]^)- 
With this formalism we will derive the complete single-spin-dependent cross section arising 
from the twist-3 three-gluon correlation functions of the nucleon. We will also clarify the 
difference of our result from that of [32] . 

The remainder of this paper is organized as follows: In section 2, we first define a 
complete set of the three-gluon correlation functions in the transversely polarized nucleon. 
We show that, actually, the three-gluon correlation functions defined in [24] are not all 
independent, i.e., some of them are redundant. So we newly define a genuine complete set of 
the three-gluon correlation functions. In section 3, we present our formalism for calculating 
the twist-3 single-spin-dependent cross section for ep^ — )■ eDX. We show that only a pole 
contribution of an internal propagator in the hard part leads to a real quantity relevant 
to the cross section, and that the tree-level Ward identities satisfied by the corresponding 
pole contributions play an essential role to give the factorized expression for the single- 
spin-dependent cross section, in terms of a complete set of the gauge-invariant correlation 
functions defined in section 2. In section 4, we present the final form of the single-spin- 
dependent cross section for ep^ — )■ eDX and discuss its characteristic features. Section 5 
is devoted to a brief summary of our result. In Appendix A, we summarize the relevant 
symmetry properties of the gluon correlation functions, which are used in section 3. 

2 Three-gluon correlation functions in the transversely 
polarized nucleon 

As a straightforward extension of the quark-gluon correlation functions discussed in, e.g., 
[26, 27], purely gluonic correlation functions can be defined as nucleon matrix elements of 
the three gluon fields on the lightcone [24] . Due to the different ways for contracting the 



color indices of the three gluon fields to obtain the color-singlet operators, one can define 
the two-types of correlation functions as 
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where F°" = F^f^n^ with F^"^ = (9"Af - df^A^ + gfabcA'^A^ being the gluon field strength 
tensor, d*'^" and f'"^"' are, respectively, the symmetric and anti-symmetric structure con- 
stants of the color SU(3) group, and we have suppressed the gauge-link operators which 
appropriately connect the field strength tensors so as to ensure the gauge invariance. p 
is the nucleon momentum, and S is the transverse spin vector of the nucleon normalized 
as S'^ = —1. We obtain the twist-3 contributions of (1) and (2), when we regard all the 
free Lorentz indices a, /3, and 7 to be transverse, and, in this twist-3 accuracy, p can be 
regarded as lightlike (p^ = 0). n is another lightlike vector satisfying p ■ n = 1, and, to be 
specific, we assume p^ = {p~^, 0, Oj_) and n'" = (0, n~ , Oj_); then, we have S'^ = (0, 0, Sj_). 

Taking into account the constraints from hermiticity, invariance under the parity and 
time-reversal transformations, and the permutation symmetry among the participating 
three gluon-fields, Ji decomposed the twist-3 contribution of (1) in terms of the real, 
Lorentz-scalar functions Oi and Oi associated with the six tensor structures as [24] 



2iM. 



N 



Oi{xi, X2)(7"^eT^"^ + Oi{x2, X2 - xi)(7^^e°^"^ + Oi{x,,xi - a;2)(7^"e^^"^ 



+Oi{xi, a;2)e"^P"^^ + Oi{x2, X2 - xi)e^^P"^° - Oi{xi,xi - a;2)e^"^"5^] , (3) 

where we have introduced the nucleon mass M^v in order to define Oi and Oi as dimen- 
sionless, and the similar decomposition of (2) was also introduced as [24] 



2iM 



N 



Ni{xi,X2)g''^e^P^'^ - N^{x2, X2 - x^^^^e"^"^ - iVi(xi, xi - X2)g'"'e^^''^ 



+iVi(xi, a;2)e"'^^"5^ - iVi(a;2, X2 - x^e^^^^^S" + N^{xi, xi - X2)e^"^"5^] , (4) 

with the other dimensionless, real functions A'^i and A^^i. In [24], the above four functions 
Oi, Oi, Ni and A^i were treated as independent twist-3 three-gluon correlation functions. 
However, the six tensor structures in (3) and (4) are not all independent and thus these 
decompositions actually define redundant correlation functions.^ To see this, we recall the 
identity, 

HU^al^pS _ na^uPpS _|_ ^i^^appS _|_ p.p^al3uS _|_ pS^ajSpu /g\ 

and contract both sides of this identity with the tensor g'^^S^PpUs. When the indices a, f3 
and 7 are associated with the transverse components, one obtains the relation. 



^This point was also noticed and briefly mentioned in [54]. 



which shows that the tensor structures associated with Oi and A''i are reexpressed by those 
associated with Oi and A^i, respectively, in (3) and (4). Accordingly, (1) and (2) should be 
decomposed into the three tensor structures. This can be formally achieved by omitting the 
terms associated with Oi and A^i in (3) and (4), respectively, and we newly define the two 
twist-3 functions 0{xi,X2) and N{xi,X2) as independent three-gluon correlation functions 
to represent (1) and (2) as 



0"^^{XI,X2] 

= 2iMN 



2iM, 



N 



0(xi,X2)(7"^e''^"^ + 0(0:2,0:2 - xO^-^e"^"^ + 0(xi,a:i - X2)(7^"e^P"^l , (7) 



N{xi,X2)g''^e^P''^ - iV(x2, X2 - x{)g^^ e'^^'"^ - N^x^xi - X2)g^''e^P''^] . (8) 



In this paper, we use these 0(xi, X2) and N{xi,X2) as constituting a genuine complete set to 
express the contribution of the three-gluon correlations to the twist-3 single-spin-dependent 
cross section for ep^ — ?■ eDX. Similarly to the decomposition (3) and (4) considered in [24], 
the independent tensor structures in (7) and (8) are parameterized by the common real func- 
tions O and A^, respectively, as consequences of the constraints from hermiticity, invariance 
under the parity and time-reversal transformations, and the permutation symmetry among 
the participating three gluon-fields. In particular, these constraints imply the following 
symmetry relations for 0{xi,X2) and N{xi,X2)i which are associated with the C-odd and 
C-even combinations of the three gluon operators in (1) and (2), respectively: 

0(xi, 0:2) = 0(0:2, xi), 0(0:1,0:2) = 0(-a:i,-X2), (9) 

N{Xi,X2) = N{x2,Xi), N{xi,X2) = -N{-Xi,-X2). (10) 

The authors of [45] also pointed out that there are only two independent three-gluon correla- 
tion functions at twist-3, in the study of the evolution equations for the twist-3 distributions 
in the transversely polarized nucleon. We also mention that the O-even three-gluon corre- 
lation function in (8) contributes to the twist-3 transverse-spin structure function g2{x, Q"^) 
in the inclusive DIS with polarized beam and target [46], and the twist-3 local operators 
associated with the double moments of (2) have been analyzed in the framework of the 
operator product expansion [47]. 

As we will see in the next section, a contribution to the single-spin-dependent cross 
section based on the corresponding factorization formula at the leading order in QCD 
perturbation theory is generated from a pole arising in the propagators in the hard partonic 
subprocesses, and such a pole contribution fixes the momentum fractions in the correlation 
functions 0"^'''(xi, X2) and A^"^'''(xi, X2) at a:i = X2 = x. This represents the situation where 
one of the three gluon-lines participating in the partonic subprocesses has zero momentum 
(see (1), (2)), and thus the corresponding contributions relevant to SSA can be referred to as 
the soft-gluon-pole (SGP) contributions. Combined with the above decompositions (7) and 



(8), the single-spin-dependent cross section proves to receive contributions associated with 
0{x, x), 0{x, 0), N{x, x) and N{x, 0). In particular, we will observe in the next section that 
the partonic hard part associated with 0{x,x) is different from the hard part associated 
with O{x,0), due to the difference in the tensor structures among the three terms in the 
right-hand side of (7). Similarly, the partonic hard part for N[x,x) is different from that 
for N{x,0), reflecting the difference in the corresponding tensor structures in (8). 

Here, we comment on the treatment of the three-gluon correlation functions adopted in 
the recent work [32] on the same phenomenon as discussed in the present paper, i.e., on the 
twist-3 mechanism to SSAs in SIDIS, ep^ — )■ eDX. The three-gluon correlation functions 
introduced by the authors of [32] read 

Ti^\^,^) = J^^^e-^^y^^gp^es,nApS\Ci^''F,^^0)Fr{y2)F^^ (11) 

with 

and it was claimed that the corresponding twist-3 single-spin-dependent cross section at 
the leading order in QCD perturbation theory can be expressed entirely in terms of these 
two functions of x (see (60) below and the discussion following this formula). However, 
Tq {x,x) in (11) can be obtained by the contraction of 0"^'^{x,x) and N"'^"'{x,x) in (1) 
and (2) with the particular tensor gjSa^S'^npi and thus can be written, using (7) and (8), as 

'^T^^\x,x) = -4Mjv {N{x,x) - iV(x,0)) , (13) 

^T^~\x,x) = -4Mn {0{x,x) + O{x,0)) . (14) 
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Using these relations, the twist-3 single-spin-dependent cross section obtained in [32] implies 
the same partonic hard parts for 0{x, x) and 0{x, 0), and similarly for N{x, x) and N{x, 0). 
Such result implied by [32] disagrees with our result as mentioned above: In the twist-3 
single-spin-dependent cross section induced by the three-gluon correlations, 0{x,x) and 
0(x, 0) are associated with the different partonic hard parts, and similarly for N{x,x) and 
A^(a;, 0), so that the corresponding cross section cannot be expressed by Tq {x,x) only. In 
this connection, we also note that (7) and (8) can be converted to give 0{x,x), O{x,0), 
N{x,x) and N{x,0) as 

<^(^' ^) = ^l^y i^9apt-ipnS - '2ga'y'^l3pns) 0°''^''{x, x), (15) 

oIVl]\! 

o 

<^(^' 0) = ^iy~ idcxpt^pns - 2gaye,3pns) 0"^''(x, x), (16) 

oik/TV 

^(^' ^) = ^IT~ i^gaiseypns - 2gaye,3pns) N'^'^^{x, x), (17) 

oIVij\[ 

1 

^(^' 0) = '^ir~ i9ai3e^pns - '2ga^e,3pns) N°''^^{x, x). (18) 

oikzjv 



From these relations, we see that the contributions of certain types of the twist-3 compo- 
nents in the three-gluon correlation functions, obtained as the contractions of 0°''^'^ {x , x) , 
N°''^'^[x,x) with the tensor ga-y^i3pns, were not taken into account in [32]. Taking into ac- 
count the components contracted with Qa-y^ppuSi as well as those contracted with gai3^ypnS, 
is required on general grounds by the Bose statistics of the gluon, and, in (15)-(18), there 
is no reason to anticipate that the latter types of components are more important than the 
former. 

3 Formalism for the twist-3 mechanism from three- 
gluon correlations 

3.1 Kinematics for ep^ — > eDX 

Here we summarize the kinematics for the SIDIS process, 

eii)+pHp,S)^e{i') + D{Ph)+X. (19) 

As noted in the last section, one can assume that the initial nucleon's momentum p is 
lightlike, p"^ = 0, in the twist-3 accuracy. But, we keep the mass rrih of the final charmed 
hadron (D-meson) as P^ = m\. The corresponding results for the case of the light-meson 
production can be obtained by the replacement ruh — )■ in all the formulae below. For the 
process (19), besides the masses of the participating particles, there are five independent 
Lorentz invariants as 

S,, = {p + ir, x,, = ^, Q^ = -q^ = -{i-e)\ Zf = ^^, qT = ^lm 

2p ■ q P ■ q 

Here, qt is the "transverse" component of q defined as 

9r = 9'+ff%f-^)p'--^P;. (21) 

\{P-Ph) P-Phj P-Ph 

satisfying qt ■ p = qt ■ Ph = 0. In the actual calculation we work in the "hadron frame" [48] 
where the virtual photon and the initial nucleon are collinear, i.e., both move along the 
z-ax.is. In this frame, specifically, their momenta q and p are given as 

9^ = (g°,9l = (0,0,0,-Q), (22) 

and, similarly, 

P^=f^,0,0,-^) , (23) 

and the outgoing D-meson is assumed to reside in the xz plane: 



The transverse momentum of the D-meson in this frame is given by PhT = ZfQx, which is 
true in any frame where the 3-momenta q and p are colhnear. The lepton momentum in 
this frame can be parameterized as 

i'^ = — (cosh-?/', sinh ■?/' cos 0, sinh ■?/' sin 0, —1) , 
i'^ = — (cosh-?/', sinh-i/'cos^, sinh-j/'sin^, 1) , (25) 






where 

cosh^ = =:::^^-l. (26) 

Q 

We parameterize the transverse spin vector of the initial nucleon S'^ as 

5^ = (0, cos $5, sin $5,0), (27) 

where $5 represents the azimuthal angle of S measured from the hadron plane. With the 
above definition, the cross section for ep^ — )■ eDX can be expressed in terms of Sep, Xbj, 
Q^, Zf, q^, (f) and $5 in the hadron frame. Note that and $5 are invariant under boosts 
in the g*-direction, so that the cross section presented below is the same in any frame where 
q and p are collinear. 

3.2 Collinear expansion and gauge invariance for the three-gluon 
contribution 

The differential cross section for ep^ — )■ eDX can be obtained as 

where Lf^^{i,i') = 2(£^£'j, + i^^'^) — Q^g^u is the leptonic tensor for the unpolarized elec- 
tron, and W'^'^{p,q,Ph) is the hadronic tensor. In the present study we are interested in 
the contribution to W^^{p, q, Ph) from the three-gluon correlation functions for the initial 
nucleon, in which c and c are created through the photon-gluon fusion process and one 
of them fragments into a D (D) meson. The fragmentation function D{z) for a c-quark 
to become the D-meson with momentum Ph is defined from the corresponding lightcone 
correlation function as 



T.lJ^^-''^'{mmDiPH)X){DiPh)X\^,i\w)\0) = (^^ + me),^Diz) + ..., (29) 



where the ellipses denote the terms associated with the gamma matrix structures which are 
irrelevant for the present purpose. In the left-hand side, we have suppressed the gauge-link 
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operators to be connected to the quark fields, as well as the trace over the color indices. 
Here, z is the relevant momentum fraction, w is the lightlike vector of 0(1/(5) satisfying 
Ph ■ w = 1, and Pc is the momentum of the c (or c) quark with mass rric, such that 
Pc — Phi ^ + ^w^) with r = (m^z — m\/ z)/2 to satisfy pi = ml. At the leading twist-2 
accuracy for the quark-fragmentation process, we set w^ = p^ /{Ph ■ p), and, in the hadron 
frame, pc is expressed as 



zQ f, , ^ , ml 2qT 1 I gT I ^l 



P? = ^ l + S + ai.^.O.-l + ^ + ^ ■ (30) 



where z = —. Note that, when we make the replacement, J2x\D{Ph)X){D{Ph)X\ — )■ 
\p'c){p'c\ in the left-hand side of (29), with \p'^) the state with a c-quark having the momentum 
p[ = Pc\^^^,, the ellipses in the right-hand side vanish and D{z) — )■ S{l/z' — l/z). The 
fragmentation function D{z) of (29) is factorized from Wf^i, as 

/dz 
—D{z)w^u{p,q,Pc), (31) 

where the summation over the c and c quark contributions is implicit. To extract the 
twist-3 effect in w^,y, one needs to analyze the diagrams of the type shown in Fig. 1. 
However, as shown in Appendix A, in the leading order with respect to the QCD coupling 
constant for the partonic hard scattering parts, the contribution of Fig. 1(a) can not give 
rise to the single-spin-dependent cross section, and only Fig. 1(b) contributes to SSA, due 
to the symmetry properties of the correlation functions of two and three gluon fields in 
the polarized nucleon. So we shall focus on the analysis of Fig. 1(b) below. The goal of 
our analysis is to show that all the contributions from Fig. 1(b) in the twist-3 accuracy 
can be expressed in terms of the gauge-invariant correlation functions 0°''^'^{xi,X2) and 
N'^^'^{xi,X2) defined as (1) and (2) in the previous section. For this purpose, we work in 
Feynman gauge and apply the collinear expansion to the hard scattering part in Fig. 1(b), 
keeping all the terms contributing in the twist-3 accuracy [27]. For simplicity in the notation, 
we shall henceforth omit the Lorentz indices n and u for the virtual photon in w^^ {p, q,Pc) 
of (31) and write it as w (p, q,Pc)- 

The contribution from Fig. 1(b) to w{p,q,Pc) can be written as 

w{p,q,Pc) = /^ I ^,Sf:,{hM,q,Pc)M^;;XhM), (32) 



(27r)4y (2 



iTX 



where 5'^j^'3v(^i5 ^2, q-iPc) is the partonic hard scattering part represented by the middle blob 
of Fig. 1(b) and M^^J^ (^^1,^2) is the corresponding nucleon matrix element (lower blob) 
defined as 

MZ\K k2) = 9 [ d'i I d'r,^'^^^^''^-'^^^\pS\AmA'Av)A':,mpS). (33) 



Note that, for later convenience, we include one QCD coupling constant in the definition 
of this nucleon matrix element. In (32), a real contribution relevant to the cross section for 
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(a) (b) 

Figure 1: Generic diagrams for the hadronic tensor of ep^ — )■ eDX induced by the gluonic 
effect in the nucleon. Each one is decomposed into the nucleon matrix element (lower blob), 
D-meson matrix element (upper blob), and the partonic hard scattering part by the virtual 
photon (middle blob). In the expansion by the number of gluon lines connecting the middle 
and lower blobs, the first two terms, (a) and (b), are relevant to the twist-3 effect induced 
by the gluons in the nucleon. 

SSA occurs from an imaginary part of the color-projected hard part C']^°'S'^l\{ki,k2,q,Pc) 
with (12), since C'^°'M'^^^ iki,k2) are pure imaginary quantities as shown in Appendix A. 
This means that only the pole contribution produced by an internal propagator in the hard 
part can give rise to SSA. 

In the leading order with respect to the QCD coupling constant, we find that four 
topologically distinct diagrams shown in Fig. 2, together with their mirror diagrams, give 
rise to the "surviving" pole contributions; here, a short bar indicates the quark propagator 
that produces the corresponding pole contribution. The other pole contributions turn out to 
cancel among themselves after summing the contributions of all the leading-order diagrams 
for (32). With the assignment of the momenta ki and k2 of gluons as shown in Fig. 2, 
the condition for those poles is given by {pc — k2 + kiY — rnl = 0. After we perform the 
collinear expansion and reach the coUinear limit, ki — )■ Xip {i = 1, 2), with xi,a;2 and 
X2 — Xi representing the longitudinal momentum fractions of the relevant three gluons, this 
condition reduces to Xi = X2 and hence a pole of such type is referred to as the soft-gluon 
pole (SGP). In the following, we assume that S'^l\ (^i? ^2, (l,Pc) in (32) represents the sum of 
the contributions of the diagrams in Fig. 2 and their mirror diagrams, in which the barred 
propagator is replaced by its pole contribution. Also, for simplicity of notation, we suppress 
the color indices a, b, c and the momenta q and Pc in the hard part S'^'^\ (/ci, k2, q, Pc), writing 
it simply as Sfj,,yx (/ci, ^2), and correspondingly, we write Mj^l^^ (/ci, /C2) as M^'^'^{ki, /C2). 

To perform the collinear expansion, we decompose the relevant gluon momenta ki {i = 
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h - h 






Figure 2: Feynman diagrams for the partonic hard part in Fig. 1(b), representing the 
photon-gluon fusion subprocesses that give rise to the "surviving" pole contribution for 
ep^ — )■ eDX in the leading order with respect to the QCD coupling constant. The short bar 
on the internal c-quark line indicates that the pole part is to be taken from that propagator. 
In the text, momenta are assigned as shown in the upper- left diagram, where Pc denotes 
the momentum of the c-quark fragmenting into the D-meson in the final state. The mirror 
diagrams also contribute. 



1,2) as 



fcf = {ki ■ n)p^ + {h ■ p)n^ + K= Xip" + oo^.kl 



(34) 



where Xi = ki ■ n and cu'*^ = g^^ — p^n,j. Since p^ ~ 5'+<5 i^ the hadron frame with (23) 
and thus the component along p^ gives the leading contribution in (34) with respect to the 
hard scale Q, we expand S^^xiki^ ^2) around ki = xip. Expressing also the gluon field A" 
in the Feynman gauge as 



A" = (p"n^ + w^J A^ =p''n-A + u;"^A^ 



(35) 



we note that, in the matrix element M'^'^^{ki, ^2) of (33), the components associated with the 
second term of (35) give rise to the contributions suppressed by ~ 1/Q or more, compared 
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with the corresponding contribution due to the first term, p°'n-A (see, e.g., [27]). According 
to the decomposition (35), the integrand of (32) can be expressed as 

S^,x{h,k2)M^''\h,k2) 

= ^ppp (A;i, k2) M"""(A;i, k2) + 5«pp {ki, k2) a;\M'^""(A;i, k2) + S.,^, {k,, k2) u\M^''^{k,, k2) 
+ ■■■ + 3,^^ (A;i, k2) u\u''^M''''^{k^, k2) + S^^, {k,, k2) u\u^^u\M''-'^{k,, k2), (36) 

where Sppp{kuk2) = S^,x{ki,k2)p^'p''p\ M"""(A;i, A;2) = M^'''\kl,k2)n^n,nx, etc. By per- 
forming the collinear expansion of the hard part of each term in the right-hand side of this 
formula, we can organize the integrand of (32) based on the order counting with (34), (35), 
keeping the terms necessary in the twist-3 accuracy. For the first term in the right-hand 
side of (36), the Taylor expansion about kt = Xip gives, 

Jppp[l^ij k2) 



Sppp{xi, X2) + uj'^f^ki ■ 



dk^ 



., ,a TK dSppp{ki,k2) 



fCi — XiP 



dk^ 



Ki — ^iP 



2 " ' " ' dkfdk^. 



/Cj — XiP 



2 " ' " ' dk?dk^ 



Ki — XiP 



+f.;" k'' I J k^ ^ Sppp{ki,k2) 
+"" ''"' "" ^^' dkfdk^2 



Ki — X-ip 



--uj\k^,uj^,kluj\kl ^'^PPP^^^^^^^ 
6 " ' " ' " ' dkfdk^.dkj 

-]-uj\kiuj^kiuj\ki ^'^^^^(^;-^^) 

2 " ' " ' " ' dkfdk^dk] 



Ki — XiP 



1 

—c 
6 



+ ^io\k!^io'',klco\k^2 



9 Sppp{ki, k2) 



dki^dk^dk^ 



I 1^" U^^t^ Ur^l L<T d^Sppp{ki,k2) 



Ki — X^P 



dkfdk^dkj 



+ 



K'i — XiP 
Ki^XiP 

(37) 



Here and below, we use the notation S^^x{xi,X2) = S^y\{xip,X2p) for simplicity. We 
have written down the expansion explicitly up to the third-order terms. This is because, 
according to the above notice with the decompositions (34), (35), the third-order terms in 
(37) behave as the same order as the first term in 

5a/3-y(^l, k2)u''^U^,u\M^^'^{ki, k2) = [S^p-^^i, X2) + ■ ■ "] u'^ ,U^ ,Uj\M^^'^ {k^, k2), (38) 

which is obtained by the Taylor expansion of the last term of (36). If we substitute (38) 
directly into the integrand of (32) and perform the integrals over ki and /c2, the first term 
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(a) (b) 

Figure 3: Ward identities used for the pole contribution. The dotted hne represents a 
scalar-polarized gluon, and the quark lines marked by a bar are on-shell. 

in the right-hand side produces the contribution, which is associated with the hard scat- 
tering between the three physical gluons from the nucleon and behaves as the same or- 
der as the formal convolution of Sai3'y{xi,X2) with the twist-3 correlation functions of (1) 
and (2). Actually, this corresponds to a quantity of twist-3. Compared to this, the el- 
lipses in (38) give rise to the terms suppressed by 1/Q or more corresponding to twist-4 
and higher, and thus are irrelevant here. We can write down the collinear expansions 
for the hard parts associated with the other terms in (36), similarly as (37): We expand 
Sapp{ki,k2), Spap{ki,k2) and Sppa{ki,k2) through the terms involving the second deriva- 
tives, and Sai3p{ki,k2), Sapp{ki^k2) and Spaji{ki^k2) through the terms involving the first 
derivatives. Thus, the collinear expansion of Sf^ux{ki, k2) in (32) produces lots of terms as 
above, and each of those terms is not gauge invariant. At first sight, it looks hopeless to 
reorganize those into a form of the convolution with only the gauge-invariant correlation 
functions 0"^'^{xi,X2) and A^"^'^(xi,a;2) of (1) and (2) used. 

However, as was the case for the pion production associated with the twist-3 quark-gluon 
correlation functions [27], great simplification occurs due to Ward identities satisfied by the 
corresponding partonic hard-scattering function, Sfj,ux{ki, k2). To derive the Ward identities, 
we note that the contribution to S^yx{ki, ^2) from each diagram in Fig. 2 all contains the 
two delta functions, 5 {{k2 + q — Pc)"^ — rnl) and 6 {{pc + ki — /c2)^ — m^), representing the 
on-shell conditions associated, respectively, with the final-state cut on the unobserved c- 
quark line and with the unpinched pole contribution of the barred propagator. Also, in each 
diagram in Fig. 2, the c-quark line fragmenting into the final-state D-meson is on-shell, see 
(30). Due to these on-shell conditions for the diagrams in Fig. 2 and the similar conditions 
for their mirror diagrams, S^^xiki, /C2) satisfies the tree-level Ward identities, 

kiS^^^xih^h) =0, 

k2S^,ux{ki,k2) = 0, 

ik2-k^)^S^,xih,k2) = 0. (39) 

Here, the last identity and the first two identities are represented diagrammatically in 
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Figs. 3(a) and 3(b), respectively. In the collinear limit, ki — )■ Xip {i = 1, 2), we have 
^ ((Pc + ki — k2)^ — ml) — )■ (l/2pc ■ 'p)^{.^i ^ ^2)1 while the other delta function associated 
with the diagrams in Fig. 2 implies X2 > xi^j (see (59) below); the similar relations hold also 
for the corresponding mirror diagrams. Thus, the collinear expansion of Ward identities of 
(39) produces a series of relations in the collinear limit: 



Spu\{Xl-,X2) = 0, 
S^p\{Xi,X2) = 0, 
{X2 — Xi)Sfj,up{xi,X2) 

dSf,px{ki,k2] 



dkf 

Saiyx{xi,X2) + Xi 



0, 



dSpux{ki,k2] 



dSpu\{ki,k2) 



Snl3\{Xi,X2) + X2 

d S^px{ki, k2 



dkf 
dS^px{ki,k2) 



dk^ 



dk^ 



fZ-i — XiP 



Ki — XiP 



0, 



Ki — XiP 



0, 



0, 



dkfdkf 



d'^Sp^x{ki,k2) 



Ki — XiP 



X2 



d'^S^px{ki,k2) 



dk?dk'o 



dk^dk^ 
dS^px{ki,k2 



Ki — XiP 



Ki — XiP 



dkf 



(40) 
(41) 
(42) 

(43) 
(44) 
(45) 
(46) 
(47) 



Ki — XiP 



and so on, where a 7^ +, /3 7^ +. These can be used to reorganize various terms obtained 
by the collinear expansion of (36), and one obtains the following results for the contribution 
in each order: 

(i) By the relations (40), (41) and (43), the first three terms in (37), and all terms arising 
in (36) up to the order in 1/Q of those three terms, vanish. 

(ii) By the relations (42), (44)- (47), the sum of the contributions of the next higher order 
in (36), which behave as the same order as the terms involving the second derivative 
in (37), vanish (see the discussion below (52)). 

(iii) The remaining contributions in (36), behaving as the same order as the first term in 
(38), eventually yield the gauge-invariant twist-3 contribution to w{p,q,pc) in (57) 
below. 

Instead of describing in detail those rather lengthy calculations in the framework of 
the standard collinear expansion, we may employ a somewhat different approach: Among 
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the relevant relations (40)- (47), in particular, (44), (45) and (47) play a role to connect 
the two terms that are generated from the different terms in the right-hand side of (36) 
through the coUinear expansion. Namely, to obtain the gauge-invariant result, we have to 
combine the contributions from different terms in (36), and, furthermore, those terms are 
associated with different numbers of derivatives. These facts suggest an approach to apply 
directly Ward identities of (39) to the second line of (36), before the collinear expansion, 
i.e., without the expansion implied by the second equality in (36), nor the Taylor expansion 
about ki = Xip. Indeed, substituting the decomposition (34) into the first two identities in 
(39), we obtain 

Spux{ki,k2) = uj''akiS^ux{ki,k2), 

Xi 

Sf,px{ki, k2) = uj^pk^Sf^uxiki, k2), (48) 

X2 

and we apply these relations to the second line of (36), yielding 
S,,x{k,,k2)M^''\k,,k2) 

= S^^x {ki, k2) —uj^a (-^"^K + h ■ ng") 

Xi 

x^u"^ (-k^2nr + k2 ■ ng^) {p^n, + uj\)M^^'^{k,, k2), (49) 

where some factors combined with the matrix element (33) in the right-hand side can be 
reexpressed as, restoring the color indices a, b and c, 

= gj <fij d'r^e^'-^e^^^'^''-^\pS\F^\Q)A'^MFrmpS), (50) 

up to the correction terms beyond the present lowest-order calculation in QCD perturbation 
theory, so that 

S'^.A (A;i, k2) M^'\h, k2) = S,,x (A;i, k2) ^^^^^^ t'Mf "(A;i, k2) + uj\Mf"{k^, k2) . 

X\X2 '- -' 

(51) 
As the next step, we perform the collinear expansion of the hard-scattering function. 
Based on the definition (50) and the decomposition (35), we note that the first and second 
terms in the parentheses in (51) correspond, respectively, to the "second" and "third" orders 
in the order counting relevant to the collinear expansion like (37). Thus, the collinear 
expansion of (51), up to the desired order, reads (see the discussion below (38)) 



S^,x{k^,k2)M^-^{k^,k2 






S^vp (a^i, 3^2) + io ^k^ - 



dki 
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fcj — XiP 






dH 



M' 



afSn I 



A ih, h) + S^yx (a^i, 3:2) w ^MY^ih, h 



fcj — X-iP J 



(52) 



First, we consider the first term in the right-hand side, which is proportional to S^^p (xi, 0:2) = 
Spivx{xi,X2)p^ and corresponds to the second-order term in the above-mentioned order 
counting. One can show that, by the direct diagrammatic calculation of S^^p (a;i,X2), the 
corresponding SGP contributions from the diagrams in Fig. 2 cancel with those from their 
mirror diagrams (i.e., S^^^p (xi,X2) arising in (42) equals zero even for xi = X2), and thus 
the first term in the parentheses (■ ■ ■) in the right-hand side of (52) vanishes. It is worth 
noting that the similar vanishing property of the SGP contributions was used also for the 
case of the pion production associated with the twist-3 quark-gluon correlation functions 
(see, e.g., [27]). On the other hand, S^^x {xi,X2)o~!^^, arising in the last term in (52), does 
not vanish, but this can be recast using the relation. 



Sf^ux (a^i, 3:2) CO ^= {xi-X2)uj ^ 



X dS^up{ki,k2) 



dk^ 



(53) 



=XiP 



which is obtained by the collinear expansion of the last Ward identity of (39). Furthermore, 
for the second term in the right-hand side of (52), we may use another relation. 



dSfj_,yp{ki, k2 



dkt 



dS^ypiki^ k2 



dkk 



(54) 



=XiP 



which can be derived for A =_L by direct inspection of the diagrams in Fig. 2 and their mirror 
diagrams. We remind that a similar relation also holds for the hard part corresponding to 
the quark-gluon correlation functions, as discussed for the case of the pion production [27]. 
Thus, (52) yields, at the twist-3 accuracy. 



Sf,ux{k,,k2)M^''\ki,k2 



1) 
^_V^^x 



dS^up{ki, k2 



3:1X2 



dk^ 



Ki — ^iP 



X [- (k^ - k^) n, + {k, - k2) ■ ng:] Mf^^k^, k2), 
where the second line gives, similarly as in (50), 

[- (A;^ - A;«) n. + {k, - k2) ■ ng^;] Mf^k,, k2) 



(55) 



= i9 I d^i I Ae^'''^e*^'^'"'''^''(P'5|F/^"(0)F;"(r7)F^""(Ob^), (56) 

up to the higher-order corrections beyond the present accuracy. Substituting these results 
into (32), we obtain the final form for the relevant twist-3 contribution to the hadronic 
tensor in SIDIS, as the factorization formula in terms of the gauge-invariant three-gluon 
correlation functions. 



w{p,q,Pc 



dx]_ r dx2_ dS^ll{ki,k2,q,Pc)p 
xi J X2 dko 



uj>^^uj^,uj%MfX(x,,^2). (57) 



fcj — Xip 
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where we have restored the color indices a,b,c as well as momentum variables q,Pc, which 
were associated with the hard-scattering function in (32), and ■Mpabd^i^^'i) denote the 
three-gluon lightcone correlation functions, obtained by integrating (56) over A;~,kj_|_ {i = 
1, 2), as 

= ^d^''0^^'{xuX2) - ^r''N^^^{x,,X2), (58) 

with 0°''^'^{xi,X2) and N°''^'^{xi,X2) in (1) and (2). As we demonstrated above in deriv- 
ing these results, all the gauge-noninvariant terms that could potentially contribute to 
w{p,q,pc) vanished or canceled among themselves, and the total twist-3 contribution to 
w{p,q,Pc), relevant to SSA, proves to be expressed solely in terms of the gauge-invariant 
three-gluon correlation functions 0{xi,X2) and N{xi,X2) defined as (7) and (8). 

When one calculates the relevant hard part, dS'^^^i^ki, k2, q, Pc)p''^ /dk!^ , arising in 

(57), one should note that the derivative with respect to k2 can hit the delta functions, 
S ((A;2 + q — PcY — 1^1) and 6 ((pc + ki — k2Y — m^), which are involved in S'^|j,'3^(A;i, k2, q,Pc)p^ 
as mentioned above (39). Such derivatives of these delta functions can be reexpressed by 
the derivatives with respect to X2, and then be treated by integration by parts, giving rise 
to the derivative of the three-gluon correlation functions of (7) and (8). After such manip- 
ulations, the former of the above delta functions, which represents the on-shell condition 
for the final-state cut on the unobserved c-quark line in the diagrams of Fig. 2, becomes 

with X = -^, and this factor appears in the final expression for our cross section based on 
(57). 

Here, we make a brief comment on the calculation presented in [32] . Using the notation 
in the present paper, the authors of [32] calculated w {p, q, Pc) with the following formula, 
in place of the right-hand side of (57): 



r dxi r 

J Xt J 



dxi f dx2 dS'p'^{ki,k2,q,Pc)pg± 



Xi J X2 9k2j_ 

where 



(^%9a^Mp^J^^{Xi,X2), (60) 

rCi^XiP 



In (60), we can make the replacement u^ — )■ g']_ , up to the irrelevant corrections of twist-4 
and higher, and, using (61) and the property S^ga/sM-pabd^i^ ^2) = 0, implied by (58) with 
(7) and (8), we see that the three-gluon correlation functions involved in (60) are indeed 
expressed by the two types of functions of (11) after evaluating the SGP at xi = X2 explicitly. 
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Clearly, (60) used in [32] leads to a result different from the result based on our complete 
formula (57). It is straightforward to see that, if the tensor structure of the three-gluon 
correlation functions A^^q],c(^i? ^2) of (58) were assumed to be given by only one structure, 
gOiP ^-ypnS ^ our formula (57) would reduce to the formula (60), up to the corrections of twist-4 
and higher. However, such assumption contradicts with the permutation symmetry required 
by the Bose statistics of the gluon, as emphasized in section 2 and represented in (7) and 



3.3 Calculation of L^^W^" 

Using the kinematic variables defined in section 3.1, the differential cross section (28) can 
be expressed as 



d^Aa a^ 



em 



dxbjdQ'^dzfdqj'd(f) QAn'^xl^SlpQ'^ 



ZfL^''{iJ')W,,{p,q,Ph), (62) 



where aem = e^/(47r) is the QED coupling constant. We restore the implicit "free" Lorentz 
indices in (57) for the virtual photon, corresponding to fj, and z/ in w^u{p,(l,Pc) of (31) 
(see the discussion above (32)), and, substituting the result into (31), we obtain the 
three-gluon-correlation contribution to Wfj,u{p,q,Ph) in (62). To calculate the contraction 
L^'^{i,i')Wf^i^{p,q,Ph) arising in (62), we introduce the following four vectors which are 
orthogonal to each other: 

T" = ^ («" + 2xt,p") , 



Y^ 




Q' 



Q 



0^ 

Z^ = -7T- (63) 



These are the extension of four basis vectors introduced for the massless case (m/j = 0) 
in [48] to the case of massive meson with TTih in the final state [32] . These vectors become 
T^ = (1, 0, 0, 0), X'^ = (0, 1, 0, 0),Y^ = (0, 0, 1, 0), Z^ = (0, 0, 0, 1) in the hadron frame. In 
the present case, we find that W'^'^ can be expanded in terms of the following six independent 
tensors [29] : 

Vr = X^X"" + Y^'Y", V^" = g^" + Z^Z", 

V^" = T^'X" + X^T^ vr = X^'X'' - Y^'Y^ 

y^u ^ j.f,yu ^ Yt'T", vr = X^'Y'' + r^X^ (64) 



where we have followed the notation in [48] but have not shown the explicit form of the 
three tensors V^^qj among nine basis tensors because these three tensors are irrelevant for 
the expansion of our W^". We also introduce the inverse tensors V^"^ for the above Vj^'^: 



pv 



V[ 



vt:" 



v^^ 



.(2T^T'' + X^X'' + r^F 



MV'^'l 









vi:" 



vr 



rpfirph 



1 



.(x^x'^-r^r'^ 



(X'^F^ + F^X'- 



(65) 



Then, one obtains 



fc=l, •••,4,8,9 


Wp.K\ 


^Q' 


= L^,Vr/g' is given by 


vAi = 1 + cosh^-?/), 


A = -2, 


.4.3 = — cos0sinh2-?/', 


Aa = cos20sinh^-?/', 


^8 = 


— sin0sinh2-?/', 



fc=l, •••,4,8,9 



(66) 



Ag = sin20sinh^-?/'. 



(67) 



By the expansion (66), the cross section for ep^ — )■ eDX consists of the five structure 
functions associated with .4.1,2, -43, A^, As and Ag, respectively, which have different de- 
pendences on the azimuthal angle 0. 

In closing this section, we summarize the prescription established for calculating the 
twist-3 single-spin-dependent cross section that is generated by the three-gluon correlation 
functions of the nucleon: The corresponding differential cross section is given by (62) with 
the expansion in the right-hand side of (66), in which Wp^j is given as (31) using our 
factorization formula (57) for Wfj,i,{p,q,pc). 

4 Result for the twist-3 cross section for ep^ — > eDX 

Using the formalism presented above, we now obtain the leading-order QCD formula for the 
single-spin-dependent cross section in the SIDIS, ep^ — ?■ eDX, generated from the twist-3 
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three-gluon correlation functions 0{xi,X2) and A^(a;i,a;2) of (7) and (8), as 
dxhjd.Q'^dzfdq^dcf) 



o^lmO^selMN /-vr 



StcxISLQ^ V 2 



E A^r ^r ^^^i-ri-iVi ' 



fc=l, ■■■,4,8,9 



X 



■ z \ O^ 



a; 






a=c,c 



rf ^, , 20(x,x)\^^-. ( d ^, ^, 2O(a;,0)\^. 
0(a;, x) ^3^^ ) A(J^ + I 3-<^(^' 0) ^^^ I ^'^ 



dx 



X 



X 



0{,x,x) ^^^ ^ O(a:,0)^^4^ 



-Aa + 



X 



X 



\( d ^^. , 2Ar(x,x) 



A*i-(A;V(..0) 



2A^(a;,0)~ 



a; 



Aa^ 



Act,, AcTj. 

X a; " 



(68) 



where the subscript k runs over 1,2,3,4,8,9 with Ak defined in (67) and Sk defined as 
Sk = sin $5 for /c = 1, 2, 3, 4 and S^ = cos $5 for A; = 8, 9. The quark-fiavor index a can, in 
principle, be c and c, with (5c = 1 and 6c = —1, so that the cross section for the l)-meson 
production ep^ — )■ eDX can be obtained by a simple replacement of the fragmentation 
function to that for the D meson, Da{z) — )■ Da{z). ag = g^/{4:7i) is the strong couphng 
constant, and Cc = 2/3 represents the electric charge of the c-quark. The lower limits of 
the integrals are given by 



Zf 



[l-Xbj)Q' 
2xbjml 



Axijjml 



2 

Xhj Qx 



:i-^6,)Q2j 



(69) 



and 






Xbj 



Xbj 



1 + — ^ 



zlq^+rn^ 



zj{l-Zf)Q^ 



1 + ^ l + \/l + S 



for zAl + Jl 



q't 



for ^, 1 + W1 + ^ <1. 



Qt 



>1, 



(70) 
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Partonic hard cross sections Acrl 



, 4) are obtained as follows: 



(Aal 



A&l 



Q6^22|_ {Q^zil - z){l -2z + 2^2 _2x + 2x^ + I2xz(l - x)(l - z)) 
+2mlQ^x{2z{l -z)+ x{l - 8z + 8z^)) - Amlx^} , 
^^^{Q'z{l-x){l-z)-mlx}, 

wx q _ 22) {Q^z{l - x){l - z) - mix} {Q^z{l - 2£)(1 - z) - 2mlx] 



^'^S - Q5(l_^)3^3 



Aai 



32gT:£ 

Q6(l-i)3i2 



{Q2^(l _ £)(1 _ ^) _ ^2^} {Q2^(i _ ^) ^ ^2| ^ 



A(T^ = Aa\ = 0, 



(71) 



fAaf 



Aal 
Aal 



Aai 
Aal 



Q4%s,2 {Q^zil -z)il-2z + 2z^ -Ax + Ax^ + 2Axz{l - x){l - z)) 
+4:mlQ^x{2z{l -z)+ x{l - 8z + 8z^)) - Sm^x^}, 
2 Aal 



2 Aal 



IGqT'i 



:{Q^z{l - 2x){l - z) - 2mlxf, 



(36(l_i)3£2' 

^3(^^(1 - 2z){Q^zil -x)il-z)- mix), 



(72) 



Aa^ 
A<j3 



A(7? 



Aal 
Aal 



Aal 



legyf^ 



3(l-i)3i2 



{Q^z{l - 2x){\ -z)- 2m2£)(Q2(i _ 6z + 6^^) - 2m2) 



c/' 



-^^^{Q'^z{\-2x){\-z)-2mlx), 

8x 



-8mlQ^xz{l - 2£)(1 - -2) + Sm^x^} 



32gTf2 



____(g2^(l - i) + m2)(g2^(l - 2i;)(l - 5) - 2m2a;), 
- Q3(i!_V,. (l - 2z)iQ^zil - 2x)(l - z) - 2m2x), 



Q- 



|gp.(g25(l-5)+m,2), 



(73) 
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^^f = Q«(i-|3,. (Q^^(1 - 4i:)(l -z)- Amlx)iQ\l -62 + 62') - 2ml), 

^^1 = -^Bw-,iQ'^i^ - 4:^)(l - ^) - ^rnlx), 

^^3 = QF(i^(l - 25) {g^5^(l - 2)\1 + 12£ - 16x2) 

+AmlQ'x2{3 - 8£)(1 - -2) - \6m\x'} , 



^^8 = q3(Sff(1 - 25)(g25(l + 2x){\ -2) + 2mlx), 
^^9 = -Qf^(g^5(l - 5)(1 + £) + m^x), 

(74) 



where 



x = ^, z = ^. (75) 

X 2; 



The single-spin-dependent cross section (68) can be decomposed into the five structure 
functions, based on the different dependences on the azimuthal angles $5 and through 
the above-mentioned exphcit forms of A^^ and 5^, as 



rf^A 



= sin $5 {T\ + J-2 cos + J-3 cos 20) 



dxhjdQ''dzfd(^d(^ 

+ cos $5 ( j; sin + J's sin 20) . (76) 

The five independent azimuthal structures of this type have been observed also in the 
twist-3 single-spin-dependent cross section for ep^ — )■ cttX, generated from the quark-gluon 
correlation functions, as presented in [29, 35, 36]. Introducing the azimuthal angles 0/i and 
05 of the hadron plane and the nucleon's spin vector S*, respectively, as measured from the 
lepton plane, they are connected to the above $5 and as $5 = 0/i — 0^, = (ph- One can 
recast (76) into the superposition of five sine modulations with these new azimuthal angles 
as 

^ ^"^ = sm{(j)h - 0s) F^i^('^'^-<^s) + sm{2(j)h - 0s) ^^^^(^'^'•-'^s) + sin 0^ F*^^^"^^ 



dxi,jdQ''dzfdqj'd(f)i 

+ sin(30/, - 0s) F^'"(30,-<^s) + sin(0^ + 0^) p^sini^,+'t>s) ^ (77) 

with the relations: 

_psin(0h-0s) _ jr psin{24>h-4>s) 



-^2 +^4 psin0s _ --^2 + >^4 



-^3 + -^5 ™,-,ii'^,_4-^„^ —-^3 + ^5 



2 

jpsin{3(j>h-'l>s) — -^ 3 ~'~ •^ 5 ^sin(0h+0s) _ ~'^ 3 ~r '^ 5 ,^^-. 
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The azimuthal-angle dependence of the single-spin-dependent cross section derived in the 
TMD approach [18] was presented in a form similar to (77), so that the decomposition (77) 
into five structure functions is convenient to make connection with the TMD approach in 
the small- gr region. 

For completeness, we consider the unpolarized cross section for the SIDIS, ep — )■ eDX, 
and list the corresponding twist-2 contribution at the leading order in QCD perturbation 
theory, which gives an extension of the study in [49, 48, 50] to the case with massive-hadron 
production in the final state. The result is 



.2 02 ^2 7 Z^ -^M ~7 I — 2^ ^aK^)^y-^)Uk 



dxbjdQ'^dzfdqfdcj) SnxhS^ Q^ 






bj'-'ep^ ^ fc=l 

-(|-(-i)(-i)-if)' <-' 

with (67), (69), and (70). This is generated from the unpolarized gluon-density distribution 
G{x) for the nucleon, 

G{x) = -i^j f^e^'^{p\Ff"iO)Fri\n)\p), (80) 

X J Ztt 

and the partonic hard cross sections are obtained as, utilizing the formalism discussed in 
section 3.3, 

<3-f = Q4(i^,)2,2 {Q'^(l -z){l-2z + 2z^ -2x + 2£2 + 12x5(1 - x){l - z)) 

+2mlQ^x{2z{l -z)+ x{l - 8z + Sz^)) - Am^^x^}, 
^^ = Q^fBmiQ'^i^ -x){l-z)- mix), (81) 

^3 = Q^?dFl(l - 25) (Q'^(l - 2£)(1 -z)- 2mlx), 
^ af = ^,^iM_(g2^(i _s:){l-z)- mlx){Q'z{l - z) + ml). 

We note that these partonic hard cross sections coincide with the corresponding results 
presented in [32], except for 6"[^. Comparing (81) for o"^ with the above result (71) for 
Aal, which represents the partonic hard cross sections associated with the derivatives, 
dO{x,x)/dx and dN{x,x)/dx, of the three-gluon correlation functions in (68), one finds 
the following relations between the partonic hard cross sections at the twist-3 level and 
those at the twist-2 level: 

The other partonic cross sections at the twist-3 level, A&l {i = 2, 3, 4) of (72)-(74), are also 
related to the partonic cross sections (81) at the twist-2 level, although, unlike (82), the 
corresponding relations cannot be manifested by direct comparison between the formulae of 
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(72)-(74) and those of (81). These remarkable relations, as well as a single "master formula" 
behind them, will be presented elsewhere [51]. We mention that the similar master formula, 
which allows us to relate the 3 — )■ 2 partonic subprocess relevant for the twist-3 level to the 
2 — )■ 2 subprocess for the twist-2 level, was derived [29] for the case of the SGP contributions 
associated with the twist-3 quark-gluon correlation functions. 

5 Summary 

In this paper, we have investigated the single spin asymmetry for the D-meson production 
in SIDIS, generated from the twist-3 three-gluon correlation functions for the nucleon. We 
first showed, correcting the previous study in [24], that there are only two independent 
three-gluon correlation functions of twist-3, 0{xi,X2) and N{xi,X2), which correspond to 
two possible ways to construct color-singlet combination composed of three active gluons. 
Then, we have formulated the method for calculating the twist-3 single-spin-dependent 
cross section generated from the three-gluon correlations. Our formulation is based on 
a systematic analysis of the relevant diagrams in the Feynman gauge and gives all the 
contribution to the cross section at the twist-3 level in the leading order in perturbative 
QCD, guaranteeing the gauge invariance of the result. As in the twist-3 mechanism for SSA 
generated from the quark-gluon correlation functions for the nucleon, the cross section in 
the present case occurs as the pole contribution of an internal propagator in the partonic 
hard-scattering subprocess, and the corresponding contribution leads to the cross section 
expressed in terms of the four types of functions of the relevant momentum fraction x: 
0{x,x), O{x,0), N{x,x) and A^(a:,0). We find that all these four-types of functions and 
their derivatives with respect to x contribute to the final form of the cross section. These 
features discussed for the SSA in ep^ — ?■ eDX also apply to the case of An in p^p — )■ hX 
ih = TT, K, D, etc.) [52]. 

These new results have also revealed that the previous studies [32, 33] missed important 
contributions. In particular, the factorization formula for the cross section used in [32, 33] 
was written down in an ad-hoc way, and, compared to our factorization formula derived 
in this paper, involved additional projection for the Lorentz structure onto a particular 
component of the hard-scattering part as well as of the three-gluon correlation functions, 
which would have lead to the result incompatible with symmetry requirements in QCD. 
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A Symmetry constraints on the gluon correlation func- 
tions 

In this Appendix, we consider the correlation functions of the gluon fields, A^{C,), in the 
nucleon. Such correlation functions of the two- or three-gluon fields, with non-lightlike 
separations between those fields, arise in the intermediate step of the analysis of the relevant 
Feynman diagrams for the SIDIS, ep^ — )■ eDX, to the twist-3 accuracy, as discussed in 
section 3.2. We discuss the symmetry properties of those correlation functions and their 
implication on SSA: We show that Fig. 1(a) does not contribute to SSA, and that Fig. 2(b) 
can give rise to SSA. 

Similarly to (32), the contribution from Fig. 1(a) to Wfj,,y{p,q,Pc) in (31) can be written 
as, 

"^'^ Sl^Jik)MltX{k,p,S), (83) 



J'^fj.u K'^J^"{2)^ 



(27r) 

where the Lorentz indices for the virtual photon in the partonic hard-scattering part S^'fj{k) 
are suppressed for simplicity, and the corresponding nucleon matrix element M^^Jk^p^S) 
is defined as a correlation function of the type mentioned above: 



M^^{k,p, S) = j d'ie^^<{pS\Al{Q)A^Si)\pS)- (84) 

Note that the color indices of the gluon fields in this formula are summed over, since only 
the color-singlet combination is relevant as a matrix element in the color-singlet hadron. 
Also, the corresponding color projection is taken for SJ^J{k) in (83). Now, invariance under 
the parity (P) and time-reversal (T) transformations implies 

M^^^{k,p,Sr = Mi;^{k,p,-S), (85) 

and, combined with the fact that the matrix element Mjt^JkjP, S) depends on the spin 
vector S'^ linearly, we find that the spin-dependent part of M^^Ak^p, S) is a pure imaginary 
quantity. For the process ep^ — ?> eDX, the leptonic tensor L^^'^ of (28) is real. Accordingly, 
an imaginary contribution from the hard part S^J{k) is necessary to give the real contribu- 
tion to the spin-dependent cross section. However, this is impossible for the leading-order 
diagrams contributing to the middle blob in Fig. 1(a). Therefore, (83) representing the 
contribution from Fig. 1(a) does not give rise to SSA. 

Next, using the similar logic as above, we consider the contribution of Fig. 1(b) to (32). 
By projecting M^^J^ {ki,k2) in (33) into the color-singlet components, we define the two 
types of the contractions of the corresponding color indices as 

Mr^(A;i, k2,p, S) ^ Ci^'^M^'ik,, k^) 

= J dij dr^e^^-^e^'^'--'-^^{pS\C'rAmA^MKmpS), (86) 
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with Cj^"" of (12). PT-invariance implies 

M^''\h,k2,p,Sr = M^''\h,h,p,-S), (87) 

which shows that the spin-dependent parts of M!^^ (/ci, k2,p, S) are pure imaginary quan- 
tities. Accordingly, only an imaginary contribution from the corresponding color-projected 
hard part C'^"'S'^l\{ki,k2,q,Pc) can give rise to SSA with (32) in ep^ — )■ eDX. For the 
diagrams shown in Fig. 1(b), the pole contribution of an internal propagator in the hard 
part can give rise to such imaginary contribution. The corresponding unpinched poles are 
allowed in Fig. 1(b), owing to the presence of an extra gluon line connecting the hard and 
soft parts, compared with Fig. 1(a). 

At this point one might recall that the A;_|_-dependent gluon distribution function [53] 
( "gluon- Sivers function") can give rise to SSA in the framework of the TMD factorization 
and this fact may contradict with the above statement that (84) does not contribute to 
SSA. However, the gluon Sivers function can become well-defined only after supplying 
an appropriate gauge-link operator between the gluon-fields, and the gauge-link operator 
actually resums the effect of the extra gluon-lines connecting the hard and soft parts, 
i.e., the gluon-Sivers function represents some effect of the contribution from Fig. 1(b). 
Therefore, there is no contradiction between the two facts. The gluon-Sivers function may 
rather represent the same effect as the three-gluon correlation functions in the region of the 
intermediate transverse momentum, similarly to the case of the quark-Sivers function and 
the twist-3 quark-gluon correlation functions as shown in [39, 40, 41]. 
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